Phase Diagram of Two-dimensional Polarized Fermi Gas With Spin-Orbit Coupling 
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We investigate the ground state of the two-dimensional polarized Fermi gas with spin-orbit cou- 
pling and construct the phase diagram at zero temperature. We find there exist phase separation 
when the binding energy is low. As the binding energy increasing, the topological nontrivial super- 
fluid phase coexist with topologically trivial superfluid phase which is topological phase separation. 
The spin-orbit coupling interaction enhance the triplet pairing and destabilize the phase separation 
against superfluid phase. 
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1. INTRODUCTION 

The topological properties have been investi- 
gated extensively in condensed matter systems 
such as topological insulators(TIs)[l, 2], topological 
superconductors(TSCs)[3-6], etc, which are described by 
topological order [7] instead of the traditional Landau 
symmetry breaking theory. In ultracold atomic system, 
the effective spin-orbit coupling(SOC) has been realized 
recently by utilizing the spatial varying laser fields[8, 9]. 
With the technique of Feshbach resonance[10, 11], the 
spin-orbit coupled ultracold atomic systems provide a 
clean platform to investigate the topological properties 
of the condensed matter system. 

The SOC significantly changes the Fermi surface and 
largely enhances the low energy density of statc[12, 13]. 
Therefor, many interesting phases and intriguing phe- 
nomena become possible. The triplet pairing and the 
transition temperature are largely enhanced[14] while the 
pair coherence lengths are suppressed by the SOC[l-^i]. In 
three dimensions, the ground state of the Fermi system 
is enriched by the SOC[16-23]. In two dimensions, the 
superfluid phase of the spin-orbit coupled Fermi gas can 
be topologically nontrivial[24-27]. Furthermore, there is 
topologically nontrivial phase separation(TPS) which is 
the coexistence of superfluid phases with different topo- 
logical order in the trapped SOC Fermi systems with pop- 
ulation imbalance[28, 29]. 

In this paper, we investigate the uniform polarized two- 
dimensional (2D) Fermi gas with SOC near a wide Fesh- 
bach resonance at zero temperature. The phase sepa- 
ration is possible for a polarized Fermi gas without the 
SOC due to the competition between the polarization 
and the pairing interaction. To map out a exact phase 
diagram, we determined the ground state by minimizing 
the thermodynamic potential of the phase separation[30]. 



In the presence of SOC, the Fermi surface is topologi- 
cally changed. The topological phase transition(TPT) 
takes place when the excitation gap is closing. Therefor, 
the topologically nontrivial superfluid phase(TSF) shows 
up in the phase diagram against the topologically trivial 
superfluid phasc(NSF). For the phase separation phase, 
the topological phase transition much more tend to take 
place in the smaller pairing gap component state, thus 
the phase separation becomes topologically nontrivial. 

This paper is organized as follows. In Sec. 2, introduc- 
ing the Hamiltonian of 2D uniform polarized Fermi gas, 
we obtain the zero temperature thermodynamic poten- 
tial by mean field theory, and then give the gap equation 
and the number equations for superfluid phase. In Sec. 3, 
we investigate the ground state by minimizing the ther- 
modynamic potential of the phase separation phase and 
map out the phase diagram in detail. A brief conclusion 
is given in Sec. 4. 



2. FORMALISM OF THE SYSTEM 

We consider the uniform 2D polarized Fermi gas with 
SOC, which is described by the Hamiltonian: 



(1) 



where Hq is the kinetic term, Hso is the spin-orbit in- 
teraction, and Hint is the s-wave interaction between the 
two fermionic species. They take 
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where ^k,(T — hk"^ /{2m) — fi^r, cl^„ic-k,a) denotes the cre- 
ation(annihilation) operators for a fermion with momen- 
tum k and spin cr = {t,i}, A is the strength of Rashba 
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spin-orbit coupling, (^k = SLig{kx + iky), g is the bare 
s-wavc interaction strength which can be renormaUzcd 
by 



E 
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By the transformation, 
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the Eq.2 becomes 
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where ak_|_(ak.±) is the creation(annihilation) operator 
for the state with hehcity (±), ^k,± = Ck ± A/c with ^k = 
£k ^ A' the chemical potentials fj, = (fi^ + /ij_)/2, 
h = {^^ — /ij^)/2, A is the pairing potential which takes 

^ = sEk < C-k,iCk,t >• 

The Hamiltonianl can be rewritten in the helicity basis 



^'k = (ak,+ , flk,- 



as: 
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with 
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/ Ck.+ e^^k/i Ae'f^ \ 

Ae-'f^ -^k,+ e~"f^h 
\ Ae*^'- e^^i-'/i -Ck,- / 

We know that the classification of above 2D BdG 
Hamiltonian, which breaks the time-reversal symmetry 
but preserves the particle- hole symmetry, is Z class [^i]. 
The topological numbers which characterize the topologi- 
cal properties of the superfiuid phases are integer. There 
is topological phase transition at the gap closing point 
h = ^ y? 4- A^. The topologically nontrivial superfiuid 



phase show up when h,> \j [i 
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2 a_ A2. 

The Hamiltonian can be diagonalized as 

1 
2 



E ^> 



k.s^k.s'^k,; 



E (ek-i?k,s) + ^,(8) 

k,s=± ^ 



ik±(Q!k,±) is the creation(annihilation) oper- 
ator for the quasiparticles with the excitation spec- 
tra i;k,± - VCk + I A 



-A2fc2 ± 2Eo, here 



^K^{(i+ I A |2)+A2fc2. 

The thermodynamical potential is $7 = — Trln[e^''^] 
with (5 = l/{kBT). At T = 0, the thermodynamical 
potential is 

A P 

(9) 
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The pairing gap should be self-consistently determined 
with chemical potential by minimizing the thermody- 
namic potential dVL/dA = and the particle number 
equations 71^ = —dfl/dfi^. They are given as 



k k,s— ± ' ^ 
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where, n = n-^ + is the total particle number and 
p = (n-f — n^jn is the polarization. In the presence 
of SOC, the Fermi surface is topologically changed and 
the triplet pairing is possible. The condensate fraction 
should include singlet and triplet contributions = no -I- 
■n\ which are given as 



no 
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3. THE PHASE DIAGRAM IN p - XUf/Ef PLANE 

There is no guarantee that the ground state of the po- 
larized Fermi gas corresponds to one of the spatially ho- 
mogeneous states. As the competition between the pop- 
ulation imbalance and the pairing interaction, the phase 
separation becomes possible. For the polarized Fermi 
gas, the stability of the phase separation against the su- 
perfiuid should be considered like the case without SOC. 
By introducing the mixing coefficient a:;(0 < x < 1) and 
ignoring the interfaces energy between the two coexisting 
phase, the thermodynamic potential of the phase separa- 
tion can be written as 



n =^ xn{Ai) + {\ - x)n{A2), 



(14) 



where, Ai{i = 1,2) is the pairing gap of the i component 
separated state. The thermodynamic potential should 
be minimized with A^ and the mixing coefficient x. The 
number equations become n^- = xncr(Ai)-|-(l— a;)nCT(A2). 
By solving the gap equations and the number equations 
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FIG. 1: The phase diagrams in p—Xkp/EF plane with binding 
energy (a) Et = O.lEp; (b) Ei, = 0.5Ef; (c) Et = O.GEf; (d) 
Eb — I.OEf- Here, Ef = k%/2m = n-n jra. The red sohd lines 
separate the TSFs from NSFs phases. The blue dash lines 
separate the phase separation from superfluid phase. The blue 
dash-dot-dot lines are the boundaries between topologically 
trivial and nontrivial phase separation. The dot lines denote 
the A/Sf = 10"'^, above which the pairing gap is lower than 
10"^ 
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FIG. 2: The thermodynamic potential H as a function of the 
pairing gap A with the binding energy E^/Ef ~ 0.6 for (a) 
XUf/Ef = 0.4, p = 0.8; (b) XkF/EF = 0.7, p = 0.8; (c) 
XkF/EF = 0.4, p = 0.3; (d) XkF/Ep = 0.7, p = 0.3. 



selfconsistently, we construct the phase diagram in p — 
Xtp/Ep plane for different binding energy. 

First, we give the phase diagrams in p — Xkp/ Ep plane 
with different binding energy in Fig.l((a) Eb = Q.lEp', 
(b) Eb = 0.5Ef; (c) Eb = 0.6Ef; (d) Eb = l.OEp). The 
phase separation, which is coexistence of distinct topo- 
logically trivial superfluid phases, show up in the absence 
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FIG. 3: The pairing gap and the chemical potential as func- 
tions of the SOC strength Xkp/Ep with p = 0.01 for (a) 
Eb — 0.5Ef and (b) Eb — I.OEf- (c) The condensate frac- 
tions of singlet and triplet contribution as functions of the 
SOC for Eb/Ef = 0.5, 1 with p = 0.01. The above two lines 
are singlet contributions while the others are triplet contribu- 
tion. The triplet contributions are enhanced by the SOC. 



of SOC. When the polarization is larger than 0.32, the 
phase separation can not sustain against topologically 
trivial superfluid (NSF) in the phase diagram without 
SOC for Eb = O.lEp case. The critical polarization in- 
crease with the binding energy as shown in Fig.l. This 
consist with the recent result without the SOC [31]. 

In the presence of SOC, the Fermi surface is topolog- 
ically changed and other interesting topologically non- 
trivial phases are possible. There is topological phase 
transition when the excitation gap closing at the critical 
point h = sj + . The topological phase transition 
tend to take place in the high polarization area in which 
the pairing gap is low and the imbalance of the chemical 
potential is large. Therefor, the phases are TSF in the 
phase diagrams with high polarization as shown in Fig.l. 

For the phase separation phase, the topological phase 
transition much more tend to take place in the low 
pairing gap component state. The phase separation 
become topologically nontrivial when the low pairing 
gap component state become topologically nontrivial as 
shown in Fig.l(b),(c),(d). As the binding energy in- 
creasing, the topological phase separation is more pos- 
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sible. The entire phase separation is topologicahy triv- 
ial with Ef, = 0.1i?F(shown in Fig. 1(a)) and nontrivial 
with Eb = 1.0i?F(shown in Fig. 1(d)). The boundary (the 
red sold line) between the TSF and NSF merge with the 
phase separation boundary (the blue dash line) as the 
binding energy increasing. Fig.l also show that the SOC 
destabilize the phase separation against superfluid phase. 
When the SOC strength increase to a critical value, the 
phase separation disappear. 

Second, we show the behavior of thermodynamic po- 
tential toward the pairing gap in different phase regions 
of the phase diagram for Eb/Ep ~ 0.6 in Fig. 2. The 
thermodynamic potential has two degenerate minimums 
in the phase separation regions as shown in Fig.2(a)(c). 
The two distinct superfluid phases can show up and co- 
exist in the phase diagram. The two coexistent states 
are all topologically trivial in Fig. 2(a). But, the the 
smaller component state is topologically nontrivial while 
the other is topologically trivial in Fig.2(c). The thermo- 
dynamic potential in the superfluid region only has one 
minimum as shown in Fig.2(b)(d). 

Finally, we show the variation of A, /i, h and the con- 
densate fractions for very low polarization [p = 0.01) 
with Eb = Q.bEp and E^ = lEp in Fig.3. The triplet 
condensate fractions are enhanced by the SOC. The SOC 
enhance the triplet pairing in virtue of the topologically 
change the Fermi surface. Therefor, the system can 
not sustain the phase separation against the superfluid 
phase as the triplet pairing increasing as well as the SOC 
strength. 

It should be point out that the gap equation diver- 
gent as the pairing gap A reduce to zero when fi < 
-(A4+4/i2)/(4A2) or ^ < min(- | h \ /2,-\^/2), hence 
there is no boundary between the normal phase and the 
superfluid phase. We map out the boundaries (dot lines) 
for A = 0.001 as shown in Fig.l. Above the curve, the 
pairing gap is A < O.OOIE'f and exponentially decreases 
as the SOC reduce to zero ["29]. 



4. CONCLUSIONS 

We construct the phase diagram for the two- 
dimensional Fermi gas with spin-orbit coupling and pop- 
ulation imbalance near a wide Fcshbach resonance. We 
map out the stability regions of the topologically triv- 
ial and nontrivial superfluid phase, and phase separation 
in detail. As the spin-orbit coupling increasing, there is 
topological phase transition. Therefor, the topologically 
nontrivial phase separation is possible. The spin-orbit 
coupling enhance the triplet pairing and suppress the 
phase separation. The phase separation can not sustain 
against superfluid phase when the spin-orbit coupling is 
large. 
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